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Abstract 

We have analyzed pion-nucleon scattering using the manifestly relativistic covariant frame- 
work of Infrared Regularization up to 0{q^) in the chiral expansion, where q is a generic 
small momentum. We describe the low-energy phase shifts with a similar quality as previ- 
ously achieved with Heavy Baryon Chiral Perturbation Theory, < 1.14 GeV. New values 
are provided for the 0{q'^) and 0{q^) low-energy constants, which are compared with pre- 
vious determinations. This is also the case for the scattering lengths and volumes. Finally, 
we have unitarized the previous amplitudes and as a result the energy range where data are 
reproduced increases significantly. 



1 Introduction 



Pion-nucleon scattering is a fundamental process involving the lightest meson and baryon. There- 
fore, at low energies it is amenable to be studied by the low-energy effective field theory of QCD, 
Chiral Perturbation Theory (CHPT) lU |2], which takes into account both the spontaneous as 
well as the explicit breaking of chiral symmetry in strong interactions. While the baryon field 
transforms linearly under the chiral group the pions do non- linearly [3]. A first step in extending 
CHPT to the systems with baryon number one was undertaken in ref. Contrary to standard 
CHPT, it was established that due to the presence of the large nucleon mass, loops do not respect 
the chiral power counting, and the lower order counterterms are renormalized because of higher 
order loops. The power counting was recovered by applying Heavy Baryon CHPT (HBCHPT) 
[5], where the heavy components of the baryon fields are integrated out [HI [7] so that manifest 
Lorentz invariance is lost. On the other hand, for some loop functions the expansion in inverse 
powers of the nucleon mass and the loop integration do not commute so that the non-relativistic 
expansion does not converge [SlITj. The recovery of the power counting, while keeping manifest 
Lorentz invariance, was achieved by the Infrared Regularization method (IR) [8j, based on the 
ideas in ref. 0. IR was extended to the multi-nucleon sector [TU] and to multi-loop diagrams [TT]. 
Another relativistic approach to baryon CHPT is the so-called extended-on-mass-shell (EOMS) 
renormalization scheme fl2[ [13]. The latter is based on removing explicitly the power counting 
breaking terms appearing in the loop integrals in dimensional regularization since they are re- 
absorbed by the finite set of low-energy counterterms up to the order the calculation is performed. 
For a recent review on baryon CHPT see ref. |14j . 

Here we focus on the application of IR CHPT methods to low-energy pion-nucleon scattering. In 
HBCHPT there is already an extensive list of detailed calculations with several degrees of precision. 
In refs. [151 El ttSj an 0{q'^) calculation is performed, with the additional inclusion of the A(1232) 
in ref. [9]. The calculation of pion-nucleon scattering was extended up to 0{q^) in ref. ^7\, while 
isospin violation (including both strong isospin breaking terms and electromagnetism) is worked 
out up to 0{q^) in ref. [IB]. The same authors also studied the influence of the A-isobar within 



the small e-expansion [H] up to 0{e^) in ref. [20j|^ Isospin breaking corrections for the pion- 
nucleon scattering lengths to 0{q^) are calculated within IR in refs. [231 121] • This is part of an 
on-going effort for providing high precision determinations of the pion-nucleon scattering lengths 
(a recent review on this issue is ref. [25]). Within IR CHPT vrA^ scattering was already considered 
in refs. [211 EZ]- Ref. [21] performed an O(g^) one-loop calculation. Its main conclusion was that 
the one-loop representation is not precise enough to allow a sufficiently accurate extrapolation 
of physical data to the Cheng-Dashen point. On the other hand, ref. [27] was interested in the 
complementary aspects of comparing IR CHPT at O(g^) to data and to previous HBCHPT studies 
[9l[T6lll7]. The conclusions were unexpected and rather pessimistic. The description obtained was 
restricted to very low center-of-mass (CM) pion kinetic energy (less than 40 MeV), such that the 
IR results badly diverge from the experimental values above that energy in several partial waves 
[27] . In comparison, the resulting phase shifts in HBCHPT [9l [16] fit pion-nucleon phase shifts 
up to significantly higher energies and then start deviating smoothly from data. Last but not 



^^Othcr chiral power-countings including the explicit A(1232) resonance are the ^-expansion [21] and the more 
recent one of ref. 1221. 
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least, ref. [27] also found an unrealistically large violation (20-30%) of the Goldberger-Treiman 
(GT) relation [28] for the pion-nucleon coupling. As we shall show, our results are somewhat 
more optimistic that those of ref. [2^ because we obtain that IR is able to describe low-energy 
pion-nucleon scattering comparable to HBCHPT at 0{q'^). Nevertheless, the caveat about the 
large violation of the GT relation in a full IR calculation at 0{q^) remains. When this calculation 
is restricted to strict 0{q^), as in HBCHPT, more realistic values around 2% are obtained for the 
violation of the GT relation. 

As a consequence of unitarity, ttN partial wave amplitudes develop a right-hand or unitarity 
cut with a branch point at the reaction threshold. The first derivative of the partial waves at 
this point is singular. Based on this, ref. [26] advocates for applying the chiral expansion to 
the subthreshold region of the vrA^ scattering amplitude where it is expected to be smoother. 
This singularity can also be avoided by applying the chiral expansion to an interaction kernel 
which, by construction, has no right-hand cut. This is the so-called Unitary CHPT (UCHPT) 
[23 EDI EB [22] ■ One of the consequences of this framework, is that the calculated vrA^ partial waves 
fulfill unitarity. We compare in this work the purely perturbative results with those obtained by 
UCHPT, with the latter being able to fit data closely up to higher energies as shown below. 
There are other methods already employed to provide unitarized nN amplitudes from the given 
chiral expansions, e.g. [331 EH EH]. We will also explore the region of the A(1232) resonance by 
including a Castillejo-Dalitz-Dyson (CDD) pole [36] in the inverse of the amplitude. The resulting 
amplitude has the same discontinuities along the right- and left- hand- cuts as the one without the 
CDD poles. Traditionally, this was a serious drawback for the bootstrap hypothesis [371 EH] since 
the source for dynamics in this approach is given precisely by the discontinuities of the amplitudes 
along the cuts. From our present knowledge based on QCD the presence of these extra solutions 
(including an arbitrary number of CDD poles) can be expected on intuitive grounds since they 
would be required in order to accommodate pre-existing resonances due to the elementary degrees 
of freedom of QCD. Alternatively, one could also include explicitly within the effective field theory 
the A(1232) resonances as a massive field [Ml SOI ttH [211 IHl 122] . 

After this introduction we give in section [2]our conventions and Lagrangians employed, including 
kinematics and equations used to project into partial waves. The calculation at the one-loop level 
up to 0{q^) is performed in section [HI where we compare directly with the expressions given in 
ref. [26] • We also present fits to the experimental data at the perturbative level and discuss the 
resulting values for the chiral counterterms, scattering lengths and volumes and the violation of 
the GT relation. The resummation of the right-hand cut by means of UCHPT is undertaken in 
section [H where we discuss the comparison with experimental data and the significant increase of 
the energy range for the reproduction of data. The conclusions are given in section [SI 

2 Prelude: Generalities, kinematics and Lagrangians 

We consider the process 7i"-{q)N(p, a; a) — j- it"-' {q')N{p' , a'; a'). Here a and a' denote the Cartesian 
coordinates in the isospin space of the initial and final pions with four-momentum q and q' , 
respectively. Regarding the nucleons, <j{<j') and a{a') correspond to the third-components of spin 
and isospin of the initial (final) states, in order. The usual Mandelstam variables are defined as 
s = {p + q)"^ = ip' + q'Y, t = {q — q'Y = {p — p'Y and u = [p — q'Y = {p' — g)^, that fulfill 
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s + t + u = 2M^ + 2m^ for on-shell scattering, with m and M^r the nucleon and pion mass, 
respectively. It is convenient to consider Lorentz- and isospin-invariant amphtudes (along our 
study exact isospin symmetry is presummed.) We then decompose the scattering amplitude as 



u{p,a). (2.1) 



Here, the Pauli matrices are indicated by Tc. In the next section we will proceed with the calcu- 
lation of A"^ and 5^ perturbatively up to 0{q^). 

In IR the Feynman diagrams for ttN scattering follow the standard chiral power counting 



l + 2L + J2yiid^ + lni-2) , (2.2) 



where L is the number of loops, Vi is the number of vertices of type i consisting of baryon fields 
(in our case rij = 0, 2) and di pion derivatives or masses. In this way, a given Feynman diagram 
for ttN scattering counts as q'^. For the calculation of pion-nucleon scattering up to 0{q'^) we 
employ the chiral Lagrangian 

CcHPT = £S + 4t) + + ^fl + , (2.3) 

where the superscript indicates the chiral order, according to eq. fl2.2p . Here, C^n^J refers to the 
purely mesonic Lagrangian without baryons and C^^^ corresponds to the one bilinear in the baryon 
fields. We follow the same notation as in ref. [26] to make easier the comparison. Then, 

4t) = (2(nX)(x+) + + • • ■ (2.4) 

where the ellipsis indicate terms that are not needed in the calculations given here. For the 
different symbols, F is the pion weak decay constant in the chiral limit and 



.2 



U 



U , u^ = iv)d^U v) , x± = M^X""^ ± MX^M . (2.5) 



The explicit chiral symmetry breaking due to the non-vanishing quark masses (in the isospin limit 
mu = md = rh) is introduced through x = 2Bo7h. The constant Bo is proportional to the quark 
condensate in the chiral limit (Olg-'g^lO) = —BqF'^6^^. In the following we employ the so-called 
sigma-parameterization where 



TTl x)2 .nix] ■ T 



Uix) = \|l-^ + ^^jr- (2.6) 
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In eq. (12 ■4p we denote by (■ ■ ■) the trace of the resulting 2x2 matrix. For the pion-nucleon 
Lagrangian we have 

£S = - -^{u,u,){^D^'D-'^p + h.c.) + - ji^^YK, u,]^ + ... , 

+ [Du, ux]m'D''D' + sym.) + h.c.) + z^{[x^,u,]D^ + h.c.) 

+ t^l^f^ (a^^([/^A, u,]u, - u,[D,, ux])D' + h.c.) 
8m 

+ ^7^5(x+)«. + ^7^5[/^.,X-] K + . . . (2.7) 

In the previous equation m is the nucleon mass in the chiral limit {rriu = rrid = 0) and the covariant 
derivative acting on the baryon fields is given by + with = [u"^ , d^u]/2. The low-energy 
counterterms (LECs) q and di are not fixed by chiral symmetry and we fit them to vrA^ scattering 
data. Again only the terms needed for the present study are shown in eq. (12. 7p . For more details 
on the definition and derivation of the different monomials we refer to refs. [161 US] • 

The free one-particle states are normalized according to the Lorentz-invariant normalization 

(p', a'- 7|p, a- 7) = 2Ep{27if5{^' - p)5...5^y , (2.8) 

where Ep is the energy of the particle with three-momentum p and 7 indicates any internal 
quantum number. A free two-particle state is normalized accordingly and it can be decomposed 
in states with well defined total spin S and total angular momentum J. For vrA^ scattering = 1/2 
and one has in the CM frame 

|7r(-p;a)iV(p,a;a)) = ^JA^Y^{ma^i\^]^J)Y['{■py\J^i^]aa) , (2.9) 

with p the unit vector of the CM nucleon three-momentum p, i the orbital angular momentum, 
m its third component and n = m + a the third-component of the total angular momentum. The 
Clebsch-Gordan coefficient is denoted by {mim2m3\jij2j3) , corresponding to the composition of 
the spins ji and (with third-components mi and m2, in order) to give the third spin j^, with 
third-component m^. The state with total angular momentum well-defined, \J^i;aa), satisfies 
the normalization condition 

(jyf ; aa'\Jfj,i; aa) = 5jj>5^>^5ut ^^^'^ da'a^a'a ■ (2-10) 

The partial wave expansion of the ttN scattering amplitude can be worked out straightforwardly 
from eq. (12. 9p . By definition, the initial baryon three-momentum p gives the positive direction of 
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the z-axis. Inserting the series of eq. (12 .Qp one has for the scattering amphtude 

{n{-p'-a')N{p',a'-a')\T\n{-p;a)N{p,a;a))=4n l^£°(z)(maV|£^ J)(Oaa|£i J)Fr(p'm£(s) 

e,m,j 

(2.11) 

where T is the T-matrix operator and Tjg is the partial wave amphtude with total angular mo- 
mentum J and orbital angular momentum i. Notice that in eq. (12. lip we made use of the fact 
that (z) is non-zero only for m = 0. Recall also that because of parity conservation partial 
wave amplitudes with different orbital angular momentum do not mix. From eq. (12. lip it is 
straightforward to isolate Tji with the result 

Tje{a', a'; a, a) = ^ — j_ — ^tTT XI / dp {T^i-p-,a)N{p ,a';a')\T\7i{-p;a)N{p,a;a)) 

X {ma'a\ih)Yripr ■ (2-12) 

In the previous expression the resulting Tj£ is of course independent of choice of a. 
The relation between the Cartesian and charge bases is given by 



TT 



1 



\A = k') • (2.13) 

According to the previous definition of states |7r+) = — 11, -|-1), |7r~) = |1, —1) and |7r°) = |7r^) = 
|1,0), where the states of the isospin basis are placed to the right of the equal sign. Notice the 
minus sign in the relationship for |vr+). Then, the amplitudes with well-defined isospin, 1 = 3/2 
or 1/2, are denoted by Tjji and can be obtained employing the appropriate linear combinations 
of Tji{a' , a'; a, a), eq. (I2.12p . in terms of standard Clebsch-Gordan coefficients. 

Due to the normalization of the states with well-defined total angular momentum, eq. (I2.10p . 
the partial waves resulting from eq. (12.121) with well defined isospin satisfy the unitarity relation 

ImT,j, = -^|T,j,|2 (2.14) 

for IpI > and below the inelastic threshold due the one-pion production at |p| ~ 210 MeV. 
Given the previous equation, the S'-matrix element with well defined /, J and denoted by Sjji, 
corresponds to 

Sjj, = l + i^^Tin, (2.15) 

satisfying Sij^S^j^ = 1 in the elastic physical region. In the same region we can then write 

SiM = e^'''^' , (2.16) 

with 5iji the corresponding phase shifts. 
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3 Perturbative calculation and its results 



From eq. (12.21) the leading order contribution to vrA^-scattering has u = 1 and it consists only 
of the lowest order pion-nucleon vertices with di = 1 and of no loops (L = 0). These diagrams 
correspond to the first two topologies shown from left to right in the first line of fig. [H where all 
the diagrams up-to-and-including 0{q^) are shown. The 0{q'^) or next-to-leading order (NLO) 
contribution still has no loops (L = 0) and contains an 0{q'^) vertex with di = 2. It is shown by 
the third diagram in the first line of fig. [H The NLO pion-nucleon vertex is depicted by the filled 
square. The 0{q^) or next-to-next-to-leading (N^LO) contributions consists of tree-level diagrams 
with at least one vertex of di = 3 type, with the other ones with di = 1. They are shown by the 
diagrams in the second line of fig. [H where the diamond corresponds to the di = 3 vertex. Finally, 
at N^LO one also has the one loop (L = 1) diagrams involving only vertices with di = 1 from 
the LO pion-nucleon Lagrangian, C^^jf, and with di = 2 from the LO pure mesonic Lagrangian, 
a^4. The one-loop diagrams are the rest of those shown in the figure and are labelled with a latin 
letter (a)-(v). In addition, one also has the wave function renormalization of pions and nucleons 
affecting the LO contribution. The calculation is finally given in terms of m, F^^ and qa-, which 

O 

implies some reshuffling of pieces once the constants m, F and g in the chiral limit are expressed 
in terms of the physical m, F^^ and qa making use of their expressions at 0{q^) |26]. In this work 
we employ the numerical values F^ = 92.4 MeV, = 139 MeV, rriN = 939 MeV, gA = 1.267 
and fi = rriN. 

The set of diagrams in fig. [T]was evaluated within IR CHPT in ref. [26] and we have re-evaluated 
it independently. We keep the same labelling for the one-loop diagrams as in this reference for 
easier comparison. We agree with all the one-loop integrals given in detail there. Regarding their 
contributions to and we also agree with all of them except for the contributions of the 

(2) 

so-called integral , that results from the tensor one-loop integrals with one meson and two 
baryon propagators, see appendix C of ref. [26]. We find that systematically all its contributions 
as given in ref. [26] should be reversed in sign. These contributions appear in diagrams (c) + (d), 
(g)-|-(h) and (i). Apart from the direct calculation, we have checked that the expressions given 
in ref. [26] violate perturbative unitarity. The latter results because unitarity, eq. fl2.14p . is a 
non-linear relation that mixes up orders in a power expansion. Denoting with a superscript the 



chiral order so that Tjji = Tjj\ + Tjj\ + Tjj\ + the unitarity relation eq. f l2.14p up to 0{q^) 



implies 




(3.1) 
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(a) (b) (c) 




(u) (v) 



Figure 1: Set of diagrams for vrA^ scattering up to-and-including 0{q'^). The 0{q) diagrams are the first 
two in the first line, from left to right. The 0{q'^) contributions correspond to the third one still in the 
first line, where the dj = 2, 77-2 = 2 vertex is indicated with a square. The rest of the diagrams from the 
second line until the bottom of the figure are 0{q'^). The di = 3, rii = 2 vertices are indicated with a 
diamond. The one-loop diagrams only have lowest order pion-nucleon vertices. 

Our expressions fulfill eq. fl3.ip. a nd those of ref. |26] also do once the sign in front of I^^ is 
reversed for all its contributionsEj Technically we follow the general procedure of refs. [8] for 
calculating within IR and we do not give any expression for the different integrals calculated here 
because they were given already in refs. [HI EHJ E] . 

Now we proceed to compare our perturbative calculation with the experimental phase shifts 

T^^The authors of ref. [35] state in page 30 that the scattering amplitude calculated obeys perturbative unitarity. 
It seems then that the difference in the sign referred above corresponds to a typo of [Mj . 
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for the low-energy data on the vrA^ S- and P-waves (which are the relevant partial waves for 
such energies.) Since our solution is perturbative one should evaluate the phase shifts in a chiral 
expansion too. From the relation between Sij£ and Tjje, eq. fl2.15p . one has 

IpI 

cos6ijisin6ije = ReTue . (3.2) 
Which implies that 6 starts at O(g^) so that up to 0{q^) one can write 

6ije = ReTiji , (3.3) 
with Tjji evaluated in the IR CHPT series (in our present case up to 0{q^).) 
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Figure 2: (Color online.) Fits to the KA85 pion-nucleon phase shifts [l5] as a function of \fs (in GeV) 
for ^/s^„^ = 1.13 GeV in IR CHPT at ©(g^) xhe KA85-1 fit corresponds to the solid curves and the 
KA85-2 fit to the dashed ones. Data points: circles are KA85 and squares WI08 data. 
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We now consider the reproduction of the irN phase shifts of the partial wave analyses of the 
Karlsruhe (KA85) group |15] and the current one of the GWU (WI08) group [16]. The fits are 
done with the full IR CHPT calculation to 0{q'^). Due to the absence of error in these analyses 
|15| H6] there is some ambiguity in the definition of the x^- Here we follow a similar strategy 
to that of ref. [31] and define an error assigned to every point as the sum in quadrature of a 
systematic plus a statistical error, 



where is the systematic error and Cr the relative one. In ref. [16] a relative error of 3% was 
taken while in ref. [31] a 5% error was considered. In the following we take for just 0.1 degrees 
and Cr = 2%. Regarding these values for the errors notice that isospin breaking corrections 
in 7iN scattering are estimated to be rather small (ref. [H] estimates for S-waves an isospin 
breaking correction < 1%.) We then consider the larger 2% value as a safer estimate for isospin 
breaking effects not taken into account in our isospin symmetric study. Notice also that the 0{q'^) 
contributions are expected to be suppressed compared with the leading term by a relative factor 
~ {M^/Af ~ (0.14/0.5)^ ~ 0.02. Although small, a finite value for helps to stabilize fits. 
Otherwise, with = 0, extra weight is given to the small energy region close to threshold, where 
the phase shifts are smaller and absolute errors decrease. Tiny differences between the calculation 
and points in the input become then exceedingly relevant. We take = 0.1 degrees since it is 
much smaller than typical values of the phase shifts and is also the typical size for the difference 
between the phase shifts of refs. [I5lll6] in the low-energy region for the Pu partial wave (compare 
the squares [15] and the circles [1^ in fig. |2l) We have also convinced ourselves that changes in 
these values for and do not affect our conclusions. 

The function to be minimized is defined in a standard way as 



with 5th the phase shift calculated theoretically. For the minimization process we employ the 
program MINUIT [18]. 

First, we discuss the reproduction of the KA85 data [15] and later the WI08 [16] ones. As a 
first strategy, we fit directly these data from threshold up to an upper value denoted by a/s^q^., 
and consider several values for y/s^^^. A data point is included every 4 MeV in -y/i. One observes 
that the P^r degree of freedom (Xd.o.f) below 1 for y/s^^^ < 1.13 GeV, and then rises fast 
with energy so that for a/s^^^ = 1.14 GeV the Xd.o.f. is 2.1 and for y/s^^^ = 1.15 GeV it becomes 
3.6. In fig. |2]we show by the solid line the result of the fit for \fs^^^ = 1.13 GeV. At the level 
of the resulting curves the differences are small when varying y/s^^^ within the range indicated 
above. A good reproduction of the data is achieved up to around ^/s < 1.14 GeV, a similar 
range of energies to that obtained in the 0{q^) HBCHPT fits of Fettes and MeiBner [16]. From 
fig. |2]one can readily see the origin of the rise in the x^ with increasing \fs^a.x- stems from 
the last points of the partial waves P33, P31 and P\\ for which the resulting curves depart from 
them, getting worse as the energy increases. The fast rising of the P33 phase shifts is due to the 
A(1232) resonance. Though the tail of this resonance is mimicked in CHPT by the LECs, its 
energy dependence is too steep to be completely accounted for at O(g^) due to the closeness of 



err 



(5) = v/^fT^ , 



(3.4) 




(3.5) 
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LEG 


KA85-1 


KA85-2 


HBCHPT 
C(g^) [16J 


HBCHPT 
Disp. [49J 


HBCHPT 


RS 

m 


Cl 


-0.71 ±0.49 


-0.79 ±0.51 


(-1.71,-1.07) 


-0.81 ± 0.12 


-1.02 ±0.06 




C2 


4.32 ± 0.27 


3.49 ± 0.25 


(3.0,3.5) 


8.43 ± 56.9 


3.32 ±0.03 


3.9 


C3 


—6.53 ± 0.33 


—5.40 ± 0.13 


(—6.3, —5.8) 


—4.70 ± 1.16 


—5.57 ± 0.05 


—5.3 


C4 


3.87 ±0.15 


3.32 ±0.13 


(3.4,3.6) 


3.40 ± 0.04 




3.7 


di + C?2 


2.48 ±0.59 


0.94 ±0.56 


(3.2,4.1) 








C?3 


-2.68 ± 1.02 


-1.10 ± 1.16 


(-4.3,-2.6) 








C?5 


2.69 ±2.20 


1.86 ±2.28 


(-1.1,0.4) 








(il4 — (^15 


-1.71 ±0.73 


1.03 ±0.71 


(-5.1,-4.3) 








C?18 


-0.26 ±0.40 


-0.07 ±0.44 


(-1.6,-0.5) 









Table 1: Columns 2-5: Values of the low-energy constants for the KA85-1 and KA85-2 fits. The Cj 
are given in GeV~^ and the di (or their combinations) in GeV~^. The renormalization scale for di{X) 
is A = 1 GeV. The interval of values obtained in [16j by fitting low-energy ttN scattering data with 
HBCHPT at 0{q^) is given in the fourth column. Other determinations are given in columns fifth |49] 
and sixth |50j . Resonance saturation estimates are collected in the last column [50j . 



the A(1232) to the ttN threshold. Indeed, this deficiency already occurred in the 0{q^) HBCHPT 
calculation of ref. [16]. However, at C(g^) the fit to data improves because of the appearance of 
new higher order LECs [T7] . 

The resulting values for the CHPT LECs are shown in the second column of table [H denoted 
by KA85-1, in units of GeV~^ and GeV^^ for the Cj and di, respectively. Note that at 0{q^) 
only the combinations of counterterms di + d2, d^, d^, du — di^ and dig appear in tiN scattering. 
The first four combinations were already explicitly shown in the expression for d^^lf, eq. fl2.7p . 
The counterterm diQ does not appear because it is re-absorbed in the physical value of the pion- 
nucleon axial- vector coupling qa-, once the lowest order g constant is fixed in terms of the former 
[26] . Under variations of a/s^^^ most of the counterterms present a rather stable behavior, with the 
C(g^) ones being the most sensitive. The change in the LECs when varying -y/i^^^ between 1.12 
to 1.15 GeV is a source of uncertainty that is added in quadrature with the statistical error from 
the fit with a/s^^^. = 1.13 GeV, which has a xio./. The central values shown correspond to 

the same fit too. We also show in the table the values obtained from other approaches at 0{q^) 
[501 [151 HQ SI], including the 0{q^) HBCHPT fit to -kN data [T6|, the dispersive analysis within 
the Mandelstam triangle of ref. [19] and the results at 0{q^) from ref. [50], that also includes an 
estimation of the 0{q^) LECs from resonance saturation (RS). Within uncertainties, our values 
for Cl, C3 and C4 are compatible with these other determinations. Instead, C2 is somewhat larger, 
which is one of the main motivations for considering other fits to vrA^ scattering following the so 
called strategy 2, as explained below. Our values are also compatible with those determined from 
the TiN parameters up to 0{q'^) in ref. [5l] that gives the intervals Ci = (—1.2, —0.9), C2 = (2.6, 4.0) 
and C3 = (-6.1,-4.4). The threshold parameters taken in this analysis are those calculated in 
ref. [15]. Regarding the 0{q^) counterterms the comparison with HBCHPT is not so clear due 
to the large uncertainties both from our side as well as from [16]. As discussed in more detail 
below, the 0{q^) contribution is typically the smallest between the different orders studied so 
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that it is harder to pin down precise values for these counterterms. Indeed, we observe from the 
second column in table [H that (is, and di^ — di^ have large errors, much larger than those of the 
0{q^) counterterms (although the error estimated for ci is also large because the fits are not very 
sensitive to this counterterm which is multiplied by the small without energy dependence.) 
Our values for the LECs c?,, again within the large uncertainties, are compatible with those of 
ref. [16]. Only di4^ — di^ is larger in our case, out of the range given in ^ISJ by around a factor 2. 

The threshold parameters for the fit KA85-1 are collected in the second column of table O 
We have evaluated the different scattering lengths and volumes by performing an effective range 
expansion (ERE) fit to our results in the low-energy region (namely, for |p| < M^(l — M'^/m?) 
which sets the range of the ERE. ll^ The error given to our threshold parameters is just statistical. 
It is so small because the values of the scattering lengths and volumes are rather stable under 
changes of a/s^^^, and LECs within their uncertainties (taking into account the correlation among 
them.) If treated in an uncorrelated way the error would be much larger. We also vary the numbers 
of terms in the ERE expansion from 3 to 5 and the slight variation in the resulting scattering 
lengths/ volumes is also taken into account in the errors given. In the last two columns of table[2l we 
give the values from the partial wave analyses of refs. |l5lll6]. Notice that the differences between 
the central values from the latter two references are larger than one standard deviation, except 
for the P33 case. The differences between the 5*31 scattering lengths and P13 scattering volumes 
are specially large. Given this situation we consider that our calculated scattering lengths and 
volumes are consistent with the values obtained in the KA85 and WI08 partial wave analyses, 
except for the P33 one for which our result is significantly larger. It is also too large compared 
with the values obtained in the 0{q^) HBCHPT fits to phase-shifts of ref. [16] . 



Partial 
Wave 


KA85-1 


KA85-2 


KA85 


WI08 




-0.100 ±0.001 


-0.103 ±0.001 


-0.100 ±0.004 


-0.084 


'^511 


0.171 ±0.001 


0.172 ±0.002 


0.175 ±0.003 


0.171 


ao_|_ 


-0.010 ±0.001 


-0.011 ±0.001 


-0.008" 


-0.0010 ±0.0012 


'^0+ 


0.090 ±0.001 


0.092 ±0.001 


0.092'^ 


0.0883 ±0.0005 




-0.052 ±0.001 


-0.051 ±0.001 


-0.044 ±0.002 


-0.038 




-0.078 ±0.001 


-0.088 ±0.001 


-0.078 ±0.002 


-0.058 




0.251 ±0.002 


0.214 ± 0002 


0.214 ±0.002 


0.194 




-0.034 ±0.001 


-0.035 ±0.001 


-0.030 ±0.002 


-0.023 



Table 2: 5-wave scattering lengths and P-wave scattering volumes in units of and M^^, respectively. 
Our results for the fits to the KA85-1 and KA85-2 are given in the second and third columns, respectively. 
The fourth column corresponds to the values of the KA85 analysis [l5]. The values for WI08 are extracted 
from ref. [46j and the errors, when given, from ref. [52] , 

" These numbers are given without errors because no errors are provided in ref. [l5]. They are deduced 
from the KA85 ones for 053^ and ag^^ ■ 



^•^Numerical problems arising for |p| — >■ prevent to calculate directly the threshold parameters as 
limp^o|p|ReT/8^Vi|p|i+2^. 
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Due to the large values for C2 and ap^^ we consider that the fit KA85-1 is not completely 
satisfactory and try a second strategy (KA85-2). As it was commented above, the rapid increase 
in the phase shifts due to the tail of the A(1232) is not well reproduced at C(g^). As a result, 
instead of fitting the P33 phase shifts as a function of energy we fit now the function tan^pgg/lpl^ 
for three points with energy less than 1.09 GeV, where Sp^^^ is the phase shifts for the P33 partial 
wave. The form of this function is, of course, dictated by the ERE and at threshold it directly 
gives the corresponding scattering volume. We take a 2% of error for these points because within 
errors this is the range of values spanned in table [2] by the KA85 and WI08 results for ap33. A 
relative error of 2% was also taken for Cr in eq. (13. 4p . The resulting values for the LECs are 
given in the third column of Table [1] and the curves for y/s^^^ = 1.13 GeV are shown in fig. [2] 
by the dashed lines, that have a Xd.o.f. — 0.86. We observe that these curves are quite similar to 
the ones previously obtained in KA85-1. Nevertheless, for the Pu partial wave the description is 
slightly worse above 1.12 GeV and it is the main contribution to the final x^- For the P33 phase 
shifts one also observes a clear difference between the two curves as the dashed line runs lower 
than the solid line. The former reproduces the standard values for the P33 scattering volume, see 
column three of Table [21 while for the latter it is larger. This is another confirmation that the 
description of the rapid rise of the P33 phase shifts at 0{q^) enforces the fit to enlarge the value 
of the resulting scattering volume. It is remarkable that now the value of the 0{q^) LEG C2 is 
smaller and perfectly compatible with the interval of values of [IE]. It is also interesting to note 
that C3 is also smaller, which is a welcome feature especially for two- and few-nucleon systems 
that are rather sensitive to large sub-leading two-pion exchange A^A^ potential that is generated 
by the inclusion of the Ci, C3 and C4 [53]. See refs. [SH ESI [56] for a thorough discussion on this 
issue for two- and few-nucleon systems. Related to this point, one has determinations of C3 and 
C4 by a partial wave analysis of the pp and np scattering data from ref. [53 with the results 



The systematic errors are not properly accounted for yet in these determinations due to the 
dependence on the matching point that distinguishes between the long-range part of the A^A^ 
potential (parameterized from GHPT) and the short-range one (with a purely phenomenological 
parameterization.) Namely, the same authors in ref. [58] considered this issue and when varying the 
matching point from 1.8 fm to 1.4 fm the LEGs changed significantly: C3 = —5.08(28) — )■ —4.99(21) 
and C4 = 4.70(70) — )■ 5.62(69) GeV^^ With respect to the C(g^) counterterms we see that the 
central values have shifted considerably compared with KA85-1. This clearly indicates that these 
LEGs cannot be properly pinned down by fitting vrA^ scattering data. Within uncertainties d^, 
(is and di^ overlap at the level of one sigma. The LEGs di + d2 and du — di^ require to take 
into account a variation of 2 sigmas. In view of this situation we consider that one should be 
conservative and give ranges of values for these latter combination of LEGs in order to make them 
compatible 



C3 = -4.78 ±0.10 GeV 
C4 = +3.96 ± 0.22 GeV 



1 



— 1 



(3.6) 



di + d2 = +0.4 . . . + 3 GeV~^ 
d^^ - di5 = -2.4 . . . + 1.75 GeV 



(3.7) 
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These values correspond to the minimum and maximum of those shown in the second and third 
columns of table [1] allowing a variation of one sigma. 

The scattering lengths and volumes for KA85-2 are collected in the third column of Table El 
They are calculated from our results similarly as explained above for the KA85-1 fit. It is remark- 
able that now the value for the P33 scattering volume is perfectly compatible with the determi- 
nations from KA85 and WI08. We see a good agreement between our 0{q^) IR CHPT results 
and the scattering lengths/volumes for KA85. Only the Pu scattering volume is slightly different, 
though the difference between the KA85 and WI08 results is significantly large for this case too. 
One also observes differences beyond the error estimated in KA85 for the P13 scattering volume 
between the KA85 and WI08 values. Ours is closer to the KA85 one. 

It is also worth emphasizing that our fits to the phase shifts of the KA85 analysis (KA85-1 and 
KA85-2), as shown in fig. [21 offer a good reproduction of the data and the worsening for higher 
energies stems in a smooth way as in 0{q^) HBCHPT [16j. This is certainly an improvement 
compared with the previous ttN study in IR CHPT to 0{q^) of ref. |23|- this latter reference, 
data could only be fitted up to around 1.12 GeV and large discrepancies above that energy, rapidly 
increasing with energy, emerged in the S'31, P13 and Pu partial waves. 

We proceed along similar lines and perform fits of type 1 and 2 to the current solution of 
the GWU group p] (WI08). These fits are denoted by WI08-1 and WI08-2, in that order. 
The resulting curves for y/s^^^ = 1.13 GeV are shown by the solid and dashed lines in fig. [3l 
respectively. One observes very similar curves to the KA85-1 and KA85-2 fits except for the Pu 
phase shifts. Here, the agreement with the WI08 data is considerably worse. This has a clear 
translation into the values of the fo^' this partial wave, which increases almost by a factor 
3, from 20 (KA85-2) to 55 (WI08-2) (the number of fitted points is 12.) It is clear from fig. El 
that IR CHPT at 0{q^) does not compare well with the Pu WI08 phase shifts even at very low 
energies, ^/s < 1.11 GeV. For the KA85 data the situation is much better, compare with fig. [2l 
Indeed, previous solutions of the GWU (and prior VPI) group had a behavior similar to that of 
KA85 for the Pu phase-shifts, e.g. the solution SMOl employed in the analysis of ref. |27] also 
using IR CHPT at 0{q^). In view of the difficulties of our study based on IR CHPT at 0{q^) 
for reproducing the Pu phase shifts of WI08 at low energies we consider advisable a revision of 
the current solution WI08 of the GWU group and the way the t]N data affect the low-energy Pu 
phase shifts in the coupled channel approach followed [59j. 

The other distinctive features when comparing strategies 1 and 2 for the fits to the WI08 data 
are similar to those already discussed for the KA85 fits. In this way, one has for WI08-1 that 
C2 and C3 have a value in modulus larger by around 1 GeV"^ than for WI08-2. Related to this, 
the P33 scattering volume is also significantly larger for WI08-1 than for WI08-2. The values of 
the fitted LECs for WI08-1 and WI08-2 are collected in the second and third columns of Table [31 
One observes that the resulting LECs at 0{q^) are quite similar between KA85-1, WI08-1, on the 
one hand, and KA85-2, WI08-2, on the other, so that within uncertainties they are compatible 
in either of the two strategies. In the third column of Table [31 we present the average of the 
LECs from our fits in Tables [H and [31 The error given for every LEG is the sum in quadrature 
of the largest of the statistical errors shown in the previous tables and the one resulting from 
the dispersion in the central values. This is a conservative procedure which recognizes that both 
strategies are acceptable for studying low-energy t^N scattering and that takes into account the 
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Figure 3: (Color online.) Fits to the WI08 pion-nucleon phase shifts [46j as a function of \fs (in GeV) 
for = 1.13 GeV in IR CHPT at 0{ct')- The WI08-1 fit corresponds to the solid curves and the 

WI08-2 fit to the dashed ones. Data points: circles are KA85 and squares WI08 data. 



dispersion in the LECs that results from changes in the data set. Within errors, the values of the 
LECs in the last column of Table [3] are compatible with those from HBCHPT at 0{(^^ ('^i4 — c^is is 
the only counterterm that differs by more than one standard deviation from the interval of values 
of ref. [16].) 

Regarding the threshold parameters we list in the second and third columns of Table H] the 
values for the S'-wave scattering lengths and P-wave scattering volumes corresponding to the fits 
WI08-1 and WI08-2 with \fs^^^ = 1.13 GeV (the same fits shown in fig. O) The procedure for 
their determination is the same as the one already discussed for KA85-1. The largest changes 
compared with the values of the KA85-1 and KA85-2 fits, respectively, occur for the a^g^ and ap^^ 
scattering length and volume, in order. The latter also shows the largest difference between the 
results of the fits following strategy 1 and 2 (a 12% of relative difference for the KA85 case and a 
9% for the WI08 one.) Nonetheless, neither of our results for ap^^, including strategy 1 and 2 KA85 
and WI08 fits, is compatible with the value of WI08 jl6], shown in the last column of Table O 
The largest difference occurs for the value of KA85-2 which is a 50% smaller than the one from 
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WI08. Coming back to the WI08 fits, we note that the isoscalar 5*- wave scattering length a^_,_ is 
now a vanishing positive number while the Pu scattering volume has decreased, and is compatible 
with the KA85 result within one sigma. We notice that the tiny errors estimated for the threshold 
parameters resulting from our fits in Tables [2] and H] are just statistical and are determined in 
the same way as explained above for the KA85-1 fit. Of course, systematic errors due to higher 
orders in the chiral expansion and different data sets taken induce larger systematic uncertainties 
than the small errors shown. In this sense, the difference between the values obtained for each 
partial wave in these columns provides a better estimation of uncertainties. We then calculate the 



average^ n of the four values for each scattering length/volume shown altogether in Tables [2] and HI 



This is given in the last column of Table |H We also see that our averaged values for the a'^j^ and 
scattering lengths are compatible with the results obtained in ref. [60], at, = 0.0015 ± 0.0022 



''0+ 



and aQ_|_ = 0.0852±0.0018 , that takes into account isospin breaking corrections in the analysis 
of recent experimental results on pionic hydrogen and pionic deuterium data 



LEG 


WI08-1 


WI08-2 


Average 


Cl 


-0.27 ±0.51 


-0.30 ± 0.48 


-0.52 ±0.60 


C2 


4.28 ±0.27 


3.55 ±0.30 


3.91 ±0.54 


C3 


-6.76 ±0.27 


-5.77 ±0.29 


-6.12 ±0.72 


C4 


4.08 ±0.13 


3.60 ±0.16 


3.72 ±0.37 


di + d2 


2.53 ±0.60 


1.16 ±0.65 


1.78 ± 1.1 


ds 


-3.65 ± 1.01 


-2.32 ± 1.04 


-2.44 ± 1.6 


d^ 


5.38 ± 2.40 


4.83 ±2.18 


3.69 ±2.93 


di4 — (ii5 


-1.17 ± 1.00 


1.27± 1.11 


-0.145 ± 1.88 


di8 


-0.86 ±0.43 


-0.72 ± 0.40 


-0.48 ±0.58 



Table 3: Fitted LECs in units of GeV-^ (a) and GeV-^ (di) for the fits WI08-1 and WI08-2 with 
V^max = 1-13 GeV. The last columns is the average of all the fits in Tables [1] and [3l 



Finally, we show in fig. H] the different chiral order contributions to the total phase shifts 
(depicted by the solid lines) for the fit KA85-1 (shown in fig. [2] by the solid lines.) The dotted 
fines correspond to the leading result, the dashed ones to NLO and the dash-dotted ones to N^LO. 
A general trend observed is the partial cancellation between the 0{q^) and 0{q^) contributions. 
For the P-waves, the cancellation is almost exact at low energies while at higher energies the 
0{q^) contribution is larger in modulus than the 0{q^) one (except for the P31 partial wave 
where the cancellation is almost exact all over the energy range shown, so that the first order 
describes well this partial wave.) For the S-waves at low energies {^/s < 1.11 GeV) the first order 
contributions dominates, though the second order one tends to increase rapidly with energy. For 
these partial waves the second order contribution is much larger than the third order one and the 
partial cancellation between these orders is weak (even both orders add with the same sign for 
5*31 at the highest energies shown.) The smallness of the third order contribution for the S'-waves 
together with the fact that it is also clearly smaller than the second order one for most of the 

^''Not the weighted average. The given errors are calculated by adding in quadrature for each LEG the largest 
of the errors in Tables [T] and [3] and the one resulting from the average of values. 
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Partial 
Wave 


WI08-1 


WI08-2 


Average 




-0.081 ±0.001 


-0.082 ±0.001 


-0.092 ±0.012 




0.165 ± 0.002 


0.167 ± 0.002 


0.169 ± 0.004 


4+ 


0.001 ±0.001 


0.001 ±0.001 


-0.005 ±0.007 


%+ 


0.082 ±0.001 


0.083 ±0.001 


0.087 ±0.005 




-0.048 ±0.001 


-0.051 ±0.001 


-0.051 ±0.002 




-0.073 ±0.001 


-0.080 ±0.001 


-0.080 ±0.006 


ap33 


0.252 ±0.002 


0.222 ±0.002 


0.232 ±0.017 




-0.032 ±0.001 


-0.035 ±0.001 


-0.034 ±0.002 



Table 4: 5- wave scattering lengths and P-wave scattering volumes in units of and respectively, 
for the fits WI08-1 and WI08-2 with y/s^^^ = 1.13 GeV. The last column corresponds to the averaged 
values of the threshold parameters of all the fits in Tables [2] and [H 



P-waves explain the difficulties to pin down precise values for the 0{q^) LECs (the (ij's), as already 
indicated above. 

The LEG dig is important as it is directly involved in the violation of the GT relation [28] . Up 
to 0{M^) one has ^ [26] 

gAm / 2M%s 
g^N = 1 

We quantify the deviation from the GT relation by 

Act = - 1 . (3.9) 

Inserting our averaged value of dig in the third column of Table [3] into eq. (13. 8p , we then find 

Agt = 0.015 ±0.018 , (3.10) 

which is compatible with the values around 2-3% that are nowadays preferred from ttN and NN 
partial wave analyses [6T| |62| 163] . In terms of the vrA^ coupling constant, from eq. (13. 8p our value 
for (iis translates in 

^^jv = 13.07 ±0.23 (3.11) 

or p = {g^NM^/Amf /vr = 0.077 ± 0.003. Within uncertainties our result at strict 0{Ml) is 
compatible at the level of one sigma with the determinations of refs. [611 |62l [63] . 

However, IR GHPT at 0{q^) gives rise to a caveat concerning the GT relation. The point is 
that the full calculation at this order (IR GHPT contains higher orders due to the 1/m relativistic 
resummation) produces a huge GT relation violation of about a 20%, similarly as in ref. |27]. For 
the evaluation of the GT relation discrepancy in our present calculations we study the 7r~p — )■ 7r~p 
scattering. We select this particular process in the charge basis of states because the crossed u- 
channel process, Ti^p — ?■ 7r+p, is purely 1 = 3/2 and thus there is no u-channel nucleon pole, which 
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requires the same quantum numbers as for the nucleon, in the isospin hmit. Otherwise the s- 
and M-channel nucleon poles overlap for some values of the scattering angle. When projecting the 
M-channel nucleon pole in a partial wave it produces a cut for m? — 2Ml + / m? < s < m^ + 2M^, 
with the branch points very close to the nucleon pole at s = m? . As a result, there is not soft way 
to calculate the residue at the s-channel nucleon pole unless the w-channel nucleon pole is removed, 
as done by considering the 7r~p — )• tt~p scattering. The latter is finally projected in the partial 
wave Pii, with the same quantum numbers as the nucleon. The ratio of the residues at the nucleon 
pole of the full 0{q^) IR CHPT partial wave and the direct (s-channel) Born term calculated with 
qa, Mj, and m at their physical values, gives us directly the ratio between the squares of the full 



pion-nucleon coupling and the one from the GT relation^J Numerically we find that the full 



calculation gives rise to a violation of the GT relation of around 20-25%, while its strict (9(M^) 
restriction is much smaller, eq. fl3.10D . Related to this one has a significant renormalization scale 
dependence on the GT violation^J In this way, for the fit KA85-1 (second column of Tabled]) at 
A = 1 GeV one has a 22% of violation of the GT relation while for A = 0.5 GeV a 15% stems. 
On the other hand, ref. [4J performed a relativistic calculation of Act directly in dimensional 
regularization within the MS — 1 renormalization scheme and obtained a natural (much smaller) 
and renormalization scale independent loop contribution to Act- It seems then that the problem 
that we find for the calculation of Aqt with IR, obtained earlier in ref. [27], is related to the 
peculiar way the chiral counting is restored in the IR approach [Ml EH] . We tentatively conclude 
that a neat advance in the field would occur once a relativistic regularization method were available 
that conserved the chiral counting in the evaluation of loops while, at least, avoided any residual 
renormalization scale dependence. 



4 Unitarized amplitudes and higher energies 

In order to resum the right-hand cut or unitarity cut we consider the unitarization method of 
refs. [m |32], to which we refer for further details. Notice that this method does not only provide 
a unitary vrA^ amplitude but also takes care of the analyticity properties associated with the right- 
hand cut. In ref. pT] this approach was used for unitarizing the 0{q'^) HBCHPT vrA^ partial waves 
from ref. [16]. However, no explicit Lorentz-invariant one- loop calculation for vrA^ scattering has 
been unitarized in the literature yet. This is an interesting point since by taking explicitly into 
account the presence of the unitarity cut the rest of the amplitude is expected to have a softer 
chiral expansion. According to ref. [32] we express T'lji as 



Tu, = ^_,] , , , (4.1) 



"i^^Note that there is no crossed Born term for n^p — > n^p and that the LO Born term in term of physical 
parameters satisfies exactly the GT relation. 
'^^Eq. (|3.10p is renormalization scale independent because the beta function for dig is zero [IB] . 
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Figure 4: (Color online.) Different chiral orders contributing to the phase shifts for the KA85-1 fit. The 
(black) dotted, (green) dashed and (blue) dash-dotted are the first, second and third order, respectively. 
The (red) solid line is the sum of all of them. 



where the unitarity pion-nucleon loop function is given by 



9is) 



(4vr) 



ai + log 



m 



Ml 



- + s m? |p 



log(s- A + 2v^|p| 



+ log(s + A + 2v^|p|) - log(-s + A + 2v^|p|) - log(-s - A + 2v^|p| 



(4.2) 



with A = Ml — m^. The interaction kernel Tiji has no right-hand cut and is determined by 
matching order by order with the perturbative chiral expansion of Tui calculated in CHPT. In 
this way, with g{s) = 0{q), one has [32] 



^(1) ^(2) ^(3) _ ^(1) ^(2) ^(3) _^(„\(rp 



(1) 
ue 



(4.3) 
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so that 



(2) _ ^(2) 

2 



rl!l = Tl% + g{s)(T;]\) , (4.4) 



and Tije = Tjjj + T}jI + Tjjl is then replaced in eq. (14. ip . Since the resulting partial wave is 
now unitary, we calculate the phase shifts directly from the relation Tjji = ^^^^e''^'-" sin5/j£ that 
follows from eqs. ( 12. 15^ and ( 12. 16^ . 

The subtraction constant ai is determined by requiring that g{s) vanishes at the nucleon mass 
s = m?. In this way the Pu partial- wave has the nucleon pole at its right position, otherwise it 
would disappear. This is due to the fact that for this partial wave TTi vanishes at s = m? so it 

2 2 

is required that g{rn?) = 0. Otherwise Ti i^, eq. ( 14. ip . would be finite at s = m^. 

Due to the closeness of the A(1232) resonance to the vrA^ threshold it is expedient to implement 
a method to take into account its presence in order to provide a higher energy description of ttN 
phase-shifts beyond the purely perturbative results discussed in section El As commented in the 
introduction we can add a CDD pole [36] in the P33 channel so as to reach the region of the 
A(1232) resonance. The addition of the CDD pole conserves the discontinuities of the partial 
wave amplitude across the cuts. A CDD pole corresponds to a zero of the partial wave-amplitude 
along the real axis and hence to a pole in the inverse of the amplitude. We then modify eq. (14. ip 
by including such a pole in T^i^^, 



Tis,= [%-,\ + ^— + g{s)] , (4.5) 
22 \ 2 2^ s — Sp I 

where 7 and sp are the residue and pole position of the CDD pole, in order, so that two new 
free parameters enter. The amplitude Tije is determined as in eq. ( 14. 4p . We also distinguish 
here between the fits to the KA85 [IS] and WI08 [IHI phase-shifts. The fits are done up to 
^/s = -\/Smax ~ l-^S GeV for all the partial waves. One cannot afford to go to higher energies 
because of an intrinsic limitation of IR CHPT. Additional unphysical cuts and poles are generated 
by the infinite order resummation of the sub- leading 1/m kinetic energy terms accomplished in 
IR [6ll ESI m]. In our case the limiting circumstance is the appearance of a pole when the 



Mandelstam variable u = OcJ When projecting in the different partial waves this singularity 
gives rise to a strong branch point at s = 2(m^ + M^) ~ 1.34^ GeV^, which indicates the onset 
of a non-physical right-hand cut that extends to infinity and that produces strong violation of 
unitarity. This translates into strong rises of the phase-shifts calculated employing eq. ( 14. 5 p for 
energies a/s > 1.26 GeV. This is why we have taken \fs^^^ = 1.25 GeV because for higher energies 
these effects are clearly visible in the calculated phase-shifts. The to be minimized is the same 



as already used for the pure perturbative study, eq. ( 13. 5p . employing also the same definition for 



'^^Many of the tensor integrals involved in the one-loop calculations of ttN scattering develop such a pole. In 
particular, it arises in the simplest scalar two-point loop function I{u), following the notation of ref. 0. 
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err (5). The resulting fits are shown in fig. [5|, where the sohd fines correspond to tfie fit of the 
KA85 data and the dashed ones to WI08. One can see a rather good agreement with data in the 
whole energy range from threshold up to 1.25 GeV, including the reproduction of the raise in the 
P33 phase shifts associated with the A(1232) resonance. The improvement is manifest in the Pu 
partial wave although some discrepancy with the WI08 data in the lower energy region remains, 
being better the agreement with KA85 phase-shifts. Compared with the perturbative treatment 
of section E] one observes a drastic increase in the range of energies for which a globally acceptable 
description of the data is achieved. 
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Figure 5: (Color online.) Fits to the KA85 and WI08 pion-nucleon phase shifts as a function of y/s (in 
GeV) employing the unitarized irN amplitudes, eq. (jl.Sp . The solid (dashed) lines correspond to the fit 
of the KA85 (WI08) data. 

The values of the resulting LECs are coHected in Table |5l We consider that the pure pertur- 
bative study of section [3] is the proper way to determine the chiral LECs. The new values in 
Table |5] do not constitute an alternative determination to those offered in Tables [1] and [3] and 
should be employed within UCHPT studies. Nonetheless, it is remarkable that the values for the 
LECs obtained are compatible with the average of values given in the fourth column of Table El 
in particular, for the 0{q'^) LECs the central values are also rather close to the fitted values in 
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Table [5l Since we have a procedure to generate the A(1232) resonance through the CDD pole in 
eq. (14 .Sp , such agreement is surprising since the contribution of this resonance to the LECs is very 
important [50] • The point is that the typical value of 7/(5 — sp) in the low-energy region studied 
in section |3]is only around a factor 2 larger in modulus than the subtraction constant ai/(47r)^ in 
eq. (14.21) . being the latter a quantity of first chiral order. As a result, at low energies, the CDD 
pole gives a contribution that can be computed as 0{q^), since the lowest order ones comes from 
— (Tjj^)^7/(s — Sp). This explains why the values of the second order LECs are preserved, despite 
having included the CDD pole. 



LEC 


Fit 


Fit 


Partial 


Fit 


Fit 




KA85 


WI08 


Wave 


KA85 


WI08 


Cl 


-0.48 ±0.51 


-0.53 ±0.48 


«S31 


-0.115 


-0.104 


C2 


4.62 ±0.27 


4.73 ±0.30 




0.152 


0.150 


C3 


-6.16 ±0.27 


-6.41 ±0.29 




-0.026 


-0.020 


C4 


3.68 ±0.13 


3.81 ±0.16 


ao_|_ 


0.089 


0.085 


di + d2 


2.55 ±0.60 


2.70 ±0.65 




-0.050 


-0.048 


d-s 


-1.61 ± 1.01 


-1.73 ± 1.04 




-0.080 


-0.075 




0.93 ± 2.40 


1.13 ±2.18 


ap33 


0.245 


0.250 


di4 — di5 


-0.46 ± 1.00 


-0.61 ± 1.11 




-0.41 


-0.039 


di8 


0.01 ±0.21 


-0.03 ±0.20 









Table 5: Fitted LECs in units GeV"^ (q) and GeV"^ (rf.) for the fits KA85 and WI08 employing the 
unitarized partial waves. We also give the scattering lengths and volumes in units of Mj^ and M~^, 
respectively. 

The values of the resulting threshold parameters with the present unitarized amplitudes are 
collected in the last two columns of Table [5l We observe that all of them are compatible with the 
averaged values given in the last column of Table HI The P33 scattering volume turns out a bit 
too high in the lines of the values obtained with the perturbative fits following strategy 1, despite 
the reproduction of the A(1232) resonance. Finally, we also mention that similarly huge values 
for the GT violation are also obtained from the unitarized amplitudes as in the pure perturbative 
treatment. Indeed, the same value for Act, eq. (13.91) . is obtained in the unitarized case for the 
same values of the LECs because g{m'^) = (there is no CDD pole in the Pu partial wave.) 

5 Summary and conclusions 

We studied elastic pion-nucleon scattering employing covariant CHPT up-to-and-including 0{q^) 
in Infrared Regularization [8j. We followed two strategies for fitting the phase shifts provided 
the partial wave analysis of refs. [l5l HH]. In one of them, instead of fitting the P33 phase-shifts, 
we considered the reproduction of the function |p|^/ tan(5p33 around the threshold region (for 
^/s < 1.09 GeV.) The rational behind this is to reduce the impact of the A(1232) when performing 
fits to data, avoiding the rapid rise of phase-shifts with energy that tends to increase the value 
of the resulting scattering volume. An accurate reproduction of pion-nucleon phase-shifts up to 
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around 1.14 GeV results. The main difference between both strategies has to do with the values of 
the 0{(f) LECs C2 and C3, that are smaller in absolute value for strategy 2 fits. As expected, the 
P33 scattering volume is also smaller for these fits and compatible with previous determinations. 
We have discussed separately the fits to data of the Karlsruhe [15] and GWU [46] groups. We 
obtain a much better reproduction of the Pn phase shifts for the former partial wave analysis. IR 
CHPT at 0{q^) is not able to reproduce the Pu phase shifts of the current solution of the GWU 
group even at very low energies. This suggests that a revision of this solution would be in order. 
The averaged values for the LECs and threshold parameters resulting from the two strategies 
and all data sets are given in the last columns of Tables H] and [3] in good agreement with other 
previous determinations. The reproduction of experimental phase-shifts is similar in quality to 
that obtained previously with 0{q^) HBCHPT showing also a smooth onset of the departure 
from experimental data for higher energies. This is an improvement compared with previous work 
P?] . In addition, we obtain a small violation of the Goldberger-Treiman relation at strict (9(M^), 
compatible with present determinations. However, the deviation from the Goldberger-Treiman 
relation is still a caveat because when all the terms in the full IR CHPT calculation at 0{q^) are 
kept the resulting discrepancy is much higher, around 20-30%. 

We have also employed the non-perturbative methods of Unitary CHPT [321 [31] to resum the 
right-hand cut of the pion-nucleon partial waves. The A(1232) resonance is incorporated in the 
approach as a Castillejo-Dalitz-Dyson pole in the inverse of the amplitude. A good reproduction 
of the phase shifts is reached for ^/s up to around 1.25 GeV. There is an intrinsic limitation in 
IR CHPT for reaching higher energies due to the presence of a branch cut at s = 2{m? + M^) ~ 
1.34^ GeV^. Above that energy strong violations of unitarity occurs due to the onset of an 
unphysical cut associated with the infinite resummation of relativistic corrections accomplished 
in IR. This also originates a strong rise of phase-shifts noticeable already for ^/s > 1.25 GeV. The 
values of the LECs at O(g^) is compatible to those obtained with the pure perturbative study. 
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